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ABSTRACT 



In this note, we consider the observational constraints on some cosmological models by using the 
307 Union type la supernovae (SNIa), the 32 calibrated Gamma- Ray Bursts (GRBs) at z > 1.4, the 
updated shift parameter R from WMAP 5-year data (WMAP5), and the distance parameter A of the 
measurement of the baryon acoustic oscillation (BAO) peak in the distribution of SDSS luminous red 
galaxies with the updated scalar spectral index n s from WMAP5. The tighter constraints obtained 
here update the ones obtained previously in the literature. 
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I. INTRODUCTION 



Recently, some observational data have been updated or have become available. In [l|, |2| , the Wilkinson 
Microwave Anisotropy Probe (WMAP) collaboration released their 5-year observational data (WMAP5). 
The data of Cosmic Microwave Background (CMB) anisotropy have been significantly improved. Also, 
in [!, Hj], the Supernova Cosmology Project (SCP) collaboration released their new dataset of type la 
supernovae (SNIa), which was called the Union compilation. The Union compilation contains 414 SNIa 
and reduces to 307 SNIa after selection cuts. This 307 SNIa Union compilation is the currently largest 
SNIa dataset. 

On the other hand, Gamma-Ray Bursts (GRBs) were proposed as a complementary probe to SNIa 
recently d, H, S @] ■ GRBs have been advocated as standard candles since several empirical GRB lumi- 
nosity relations were proposed as distance indicators. However, there is the so-called circularity problem 
in the direct use of GRBs to probe cosmology [f| ■ Recently, a new idea to calibrate GRBs in a completely 
cosmology independent manner has been proposed [13] , and the circularity problem can be solved. The 
main idea is that of the cosmic distance ladder. Similar to the case of calibrating SNIa as the secondary 
standard candles by using Cepheid variables, which are primary standard candles, we can also calibrate 
GRBs as standard candles with a large amount of SNIa. Following the calibration method proposed 
in [l(|, the distance moduli /i of 32 calibrated GRBs at redshift z > 1.4 are derived in [ll[. Now, one 
can use them to constrain cosmological models without circularity problem. See [13, HH for details. As 
argued in [1, [32|], the observations at z > 1.7 are fairly important to distinguish cosmological models 
and break the degeneracies between the parameters. In this note, we try to combine GRBs with the 
conventional datasets to constrain cosmological models. Although the number of GRBs is small and the 
systematic and statistical errors are very large so that their contribution to the constraints would be not 
so significant, this is still a beneficial exploration. 

Here, we consider the observational constraints on some cosmological models by using the 307 Union 
SNIa compiled in Q, the 32 calibrated GRBs at z > 1.4 compiled in Table I of .11], the updated shift 
parameter R from WMAP5 [l[, and the distance parameter A of the measurement of the baryon acoustic 
oscillation (BAO) peak in the distribution of SDSS luminous red galaxies [HJ with the updated scalar 
spectral index n s from WMAP5 [l|. 

We perform a \ 2 analysis to obtain the constraints on the parameters of cosmological models. The 
data points of the 307 Union SNIa compiled in [3( and the 32 calibrated GRBs at z > 1.4 compiled in 
Table I of [ll[ are given in terms of the distance modulus /i bs{zi)- On the other hand, the theoretical 
distance modulus is defined as 



where /ip = 42.38 — 5 log 10 h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 



in which E = H / Hq and H is the Hubble parameter; p denotes the model parameters. The % 2 from the 
307 Union SNIa and the 32 calibrated GRBs at z > 1.4 are given by 



where a is the corresponding la error. The parameter /zo is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over fiQ. However, there is an alternative 
way. Following [l3|, 03] , the minimization with respect to /iq can be made by expanding the xt of Eq. 
with respect to /xq as 



Hth{zi) = 5 log 10 D L (zi) + 



(1) 




(2) 




(3) 



4(p) = i 



2ijl q B + /i 2 C, 
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where 
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B(P) = £ 
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Eq. Q has a minimum for /^o = B/C at 



X»=i(p) 



B(p) 2 



(5) 
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Since x 2 min = x 2 TOjn obviously, we can instead minimize xfi which is independent of hq. Note that the 
above summations are over the 307 Union SNIa compiled in [H and the 32 calibrated GRBs at z > 1.4 
compiled in Table I of [11]. On the other hand, the shift parameter R is defined by [lol [HI] 



R = r JL 

R -^°J E(z)' 



where the redshift of recombination z* = 1090 which has been updated in [l|, and fi m o is the present 
fractional energy density of pressureless matter. The shift parameter R relates the angular diameter 
distance to the last scattering surface, the comoving size of the sound horizon at z* and the angular scale 
of the first acoustic peak in the CMB power spectrum of the temperature fluctuations [HI, [l6[ • The value 
of R has been updated to 1.710 ± 0.019 from WMAP5 [1]. The distance parameter A is given by 



where z b = 0.35. In [17], the value of A has been determined to be 0.469 (n s /0.98) -0 - 35 ± 0.017. Here the 
scalar spectral index n s is taken to be 0.960 which has been updated from WMAP5 So, the total x 2 
is given by 



where x 2 is given in Eq. (JSJ, Xcmb = ( R ~ Robs) 2 /<? 2 R and Xbao = ~ A obs ) 2 /a\. The best-fit model 
parameters are determined by minimizing the total x 2 ■ As in [18| , the 68% confidence level is determined 
by Ax 2 = x 2 — Xmin — 1A 2.3 and 3.53 for n p = 1, 2 and 3, respectively, where n p is the number of free 
model parameters. Similarly, the 95% confidence level is determined by Ax 2 = X 2 ~ Xmin — 4-0, 6.17 
and 8.02 for n p = 1, 2 and 3, respectively. 

In sections ITT1 Mil and IIV( we consider the joint constraints on single-parameter, two-parameter and 
three-parameter cosmological models respectively, by using the 307 Union SNIa compiled in [3J], the 32 
calibrated GRBs at z > 1.4 compiled in Table I of 11], the updated shift parameter R from WMAP5 [l|, 
and the distance parameter A of the measurement of BAO peak in the distribution of SDSS luminous 
red galaxies [12] with the updated scalar spectral index n s from WMAP5 [H . Note that we also present 
the constraints without GRBs for comparison. A brief summary is given in section fVl 



In this section, we consider the constraints on three single-parameter models. They are the flat ACDM 
model, the flat DGP model and the new agegraphic dark energy (NADE) model. 




(7) 
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II. SINGLE-PARAMETER MODELS 



A. 



Flat ACDM model 



As is well known, for the spatially flat ACDM model, 

E( Z ) = y/Qmo{l + Z) 3 + (1 - Ono) • 



(9) 
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It is easy to obtain the total \ 2 as a function of the single model parameter O m o- In Fig.[TJ we present the 
corresponding likelihood C oc e~ x I 2 . The best fit has Xmin = 325.522, whereas the best-fit parameter is 

fi m0 = 0.2714t°;gi[52 ( witn 1(7 uncertainty) or Q m0 = 0.2714tao297 ( witn 2cr uncertainty). 

For comparison, we also present the likelihood without GRBs in Fig. [IJ whereas the best-fit parameter 
reads O ro0 = 0.2698±{];2i§2 ( with ltT uncertainty) or Q m0 = 0.2698to;o297 ( with 2(7 uncertainty). 
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FIG. 1: The likelihood £ oc e" x2/2 for the flat ACDM model and the flat DGP model. The results for the cases 
with and without GRBs are indicated by the black solid lines and the red dashed lines, respectively. 

B. Flat DGP model 



One of the leading modified gravity models is the so-called Dvali-Gabadadze-Porrati (DGP) braneworld 
model [H, H3| , which entails altering the Einstein-Hilbert action by a term arising from large extra 
dimensions. For a list of references on the DGP model, see e.g. [U, [22j and references therein. 

As is well known, for the spatially flat DGP model (here we only consider the self-accelerating branch), 
E(z) is given by Hfl, H S3 



(10) 



(11) 



E( Z ) = v/fWl + Zf + tt rc + 

where Q rc is constant. E(z = 0) = 1 requires 

Therefore, the flat DGP model has only one independent model parameter. Notice that < f2 rc < 1/4 
is required by < f2„ l0 < 1. 

It is easy to obtain the total % 2 as a function of the single model parameter fi rc . Also in Fig. [U we plot 
the corresponding likelihood C oc e~ x I 2 . The best fit has Xmin = 345.56, whereas the best-fit parameter 
is n ra = 0.1286±2;oo57 ( with Iff uncertainty) or fl ra = 0.1286t^"g (with 2a uncertainty). From Eq. lfTT]l . 
fi TO can be derived correspondingly. 

For comparison, we also present the likelihood without GRBs in Fig. [lj whereas the best-fit parameter 
reads D, rc = 0.1289lo.oo57 ( with lcr uncertainty) or Q rc = 0.1289io;°"e (with 2cr uncertainty). 
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C. New agegraphic dark energy model 



In [23, l2J| , the so-called "new agegraphic dark energy" (NADE) model has been proposed recently, 
based on the Karolyhazy uncertainty relation which arises from quantum mechanics together with general 
relativity. In this model, the energy density of NADE is given by [H, [24| 

3n 2 m 2 , % 

P* = —^> (12) 

where m p is the reduced Planck mass; n is a constant of order unity; 77 is the conformal time 

i-f"-l4m- <13) 



a 



in which a = (1 + z) 1 is the scale factor. Obviously, f] = 1/a, where a dot denotes the derivative with 
respect to cosmic time t. The corresponding fractional energy density of NADE reads 



p q 



^ ~ 3m 2 p H 2 # 2?? 2 ■ ( 14 ) 

From the Friedmann equation H 2 = (p m + p q ) / (3m?) , the energy conservation equation p m +3H p m = 0, 
and Eqs. (fT2"]) — l[T4"]). we find that the equation of motion for Q q is given by [23LI24] 



3(1 + Z)- 1 - \yfil q 



(15) 



From the energy conservation equation p q + 3H(p q +p q ) = 0, and Eqs. (fT2"|) — (fT4"|) . it is easy to find that 
the equation-of-state parameter (EoS) of NADE is given by [H, [HJ 

Wq^ — = -l + —- ' ( 16 ) 

p q in a 

When a — > 00, Q, q — *■ 1, thus w q —1 in the late time. When a — > 0, il q — > 0, so 0/0 appears in w q 
and hence we cannot obtain w q from Eq. (|16p directly. Let us consider the matter-dominated epoch, 
i? 2 oc p m oc a^ 3 . Thus, a}/ 2 da oc dt = adr\. Therefore, r\ oc a 1 / 2 . From Eq. (|12p . p q oc a -1 . From the 
energy conservation equation p q + 3Hp q (l + w q ) — 0, we obtain w q = —2/3 in the matter-dominated 
epoch. Since p m oc a" 3 and p 9 oc a -1 , it is expected that fl q oc a 2 . Comparing w q = —2/3 with Eq. (|T6|) . 
we find that J7 g = n 2 a 2 /4 in the matter-dominated epoch as expected. For a <c 1, provided that n is of 
order unity, !1 ? C 1 naturally follows. There are many interesting features in the NADE model and we 
refer to the original papers [HI, for more details. 

At first glance, one might consider that NADE is a two-parameter model. However, as shown in (23l |. 
NADE is a single-parameter model in practice, thanks to its special analytic feature fl q = n 2 et 2 /4 = 
n 2 (l + z)~ 2 /4 in the matter-dominated epoch, as mentioned above. If n is given, we can obtain Q q (z) 
from Eq. (|15| with the initial condition fl q (zi n i) — n 2 (l + Zi„j)~ 2 /4 at any Zi n % which is deep enough into 
the matter-dominated epoch (we choose Zi n i = 2000 as in [23|]), instead of £l q {z = 0) = fl q o = 1 — Q m o 
at z = 0. Then, all other physical quantities, such as Q m (z) = 1 — £l q (z) and w q (z) in Eq. (|16p . can 
be obtained correspondingly. So, Cl m o — £l m {z = 0), fi g o = ^g(z = 0) and u^o = w 9 (z = 0) are 
not independent model parameters. The only model parameter is n. Therefore, the NADE model is a 
single-parameter model in practice. To our knowledge, it is the third single-parameter cosmological model 
besides the flat ACDM model and the flat DGP model. 

From the Friedmann equation H 2 = (p m + p q ) j (3m 2 ), we have 



E[z) 



n m0 (i + z) 3 
i-n q {z) 



1/2 

(17) 



If the single model parameter n is given, we can obtain Q q (z) from Eq. (|15[) . Thus, we get f2 m o = 
1 — rt q (z = 0). So, E(z) is at hand. Therefore, we can find the corresponding total x 2 in Eq. JHJ. In Fig.0 



G 



we plot the corresponding likelihood C oc e~ x "/ 2 as a function of n. The best fit has Xmin = 336.061, 
whereas the best- fit parameter is n = 2.802. We present the best-fit value of n and the corresponding 
derived Cl m o, Q q o and w q o with la and 2a uncertainties in Table HI Obviously, these constraints on the 
NADE model are tighter than the ones obtained in [231 ] . 

For the case without GRBs, the best-fit parameter is n = 2.808. For comparison, we also present the 
best-fit value of n and the corresponding derived Q TO 0i Q q o and w q o with la and 2a uncertainties in 
Table D 



Uncertainty 


n 


f2 m o 




WqO 


la (with GRBs) 
2cr (with GRBs) 


9 ono+ 092 
2.802« : 1 1 8 ^ 


n 27q +0.016 

97Q+ 033 


n 79-, +0.015 
u - ' zi -0.016 

n 791 +0.030 

u - ' zl -0.033 


-o.798±S;g8S 

n 7Q8+0 0° 9 
U. f »O_ .009 


la (without GRBs) 
2a (without GRBs) 


9 ona+0 092 
2.8U»_ 090 

9 one+0-186 
2.8U8_ 179 


0.278l° ;° i« 

n 978+0-033 
U.ZfO_ 030 


n 799+O.OI5 
U. I .016 

n 7 n 9 +0.030 
U. rzz_ .033 


n 7Q8+0' 004 
U. f »O_ 004 

n 7Q8+ 009 
-u.cy»_ .oo9 



TABLE I: The best-fit value of n and the corresponding derived fi m o, Q q o and w q o with lcr and 2<r uncertainties 
for the NADE model. See text for details. 




FIG. 2: The likelihood C oc e" x /2 for the NADE model. The results for the cases with and without GRBs are 
indicated by the black solid line and the red dashed line, respectively. 



III. TWO-PARAMETER MODEL 



Here, we consider the XCDM model which is a two-parameter model. In the spatially flat universe 
which contains pressureless matter and dark energy whose EoS is a constant w, the corresponding E(z) 
is given by 

E(z) = ^fi m o(l + z) 3 + (l-fi m0 )(l + z) 3 ( 1 +»'>. (18) 
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By minimizing the corresponding total \ 2 in Eq. ([5]), we find the best-fit parameters fl m o — 0.271 and 
w = —0.951, while Xmin — 324.821. In Fig.O we present the corresponding 68% and 95% confidence level 
contours in the Q m o — w parameter space for the XCDM model. For comparison, the best-fit parameters 
are tt m o = 0.270 and w = —0.954 for the case without GRBs, whereas the corresponding confidence level 
contours are also presented in Fig. [3] 




FIG. 3: The 68% and 95% confidence level contours in the Q m o — w parameter space for the XCDM model. 
The results for the cases with and without GRBs are indicated by the black solid lines and the red dashed lines, 
respectively. The best-fit parameters for the cases with and without GRBs are also indicated by a black solid 
point and a red box, respectively. 

IV. THREE-PARAMETER MODEL 

Now, we consider the familiar Chevallier-Polarski-Linder (CPL) model [25|, (26)], in which the EoS of 
dark energy is parameterized as 

Z 

Wde = Wq + W a (l - a) = Wq + W a — , (19) 

1 + Z 



where wq and w a are constants. As is well known, the corresponding E(z) is given by [18L I27L [2 

3w a z 



E(z) 



fimoU + z f + (! - «mo) (1 + z) 3 V +w ° +w ^ exp 



1 + z 



1/2 

(20) 



There are 3 independent parameters in this model. By minimizing the corresponding total \ 2 in Eq. (jHJ , 
we find the best-fit parameters O m0 = 0.280, w — —1.146 and w a — 0.894, while Xmin = 322.475. In 
Fig. [H we present the corresponding 68% and 95% confidence level contours in the wo — w a plane for the 
CPL model. Also, the 68% and 95% confidence level contours in the f2 m o — plane and the fl m o — w a 
plane for the CPL model are shown in Fig. El It is easy to see that these constraints on the CPL model 
are much tighter than the ones obtained in |29j |. 

For comparison, the best-fit parameters are fl m o = 0.278, wq = —1.140 and w a — 0.859 for the case 
without GRBs, whereas the corresponding confidence level contours are also presented in Figs. [Hand El 
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FIG. 4: The 68% and 95% confidence level contours in the wq — w a plane for the CPL model. The results for the 
cases with and without GRBs are indicated by the black solid lines and the red dashed lines, respectively. The 
best-fit parameters for the cases with and without GRBs are also indicated by a black solid point and a red box, 
respectively. 

V. SUMMARY 

In this note, we consider the observational constraints on some cosmological models by using the 307 
Union SNIa compiled in the 32 calibrated GRBs at z > 1.4 compiled in Table I of [ll|, the updated 
shift parameter R from WMAP5 and the distance parameter A of the measurement of the baryon 
acoustic oscillation (BAO) peak in the distribution of SDSS luminous red galaxies with the updated 
scalar spectral index n s from WMAP5 w. T he tighter constraints obtained here update the ones obtained 
previously in the literature (e.g. [H, [Tj/IH, [H, [H, HU). 

It is worth noting that GRBs are potential tools which might be powerful to probe the cosmic expansion 
history up to z > 6 or even higher redshift. Of course, due to the large scatter and the lack of a large 
amount of well observed GRBs, there is a long way to go in using GRBs extensively and reliably to probe 
cosmology. The cosmology independent calibration method of GRBs proposed in [lfj| is an important step 
towards this end. The works of this note and are beneficial explorations. Along with the accumulation 
of well observed GRBs with much smaller errors, we believe that a bright future of GRB cosmology is 
awaiting us. Combining the calibrated GRBs with other probes such as SNIa, CMB, large-scale structure 
and weak lensing, we can learn more about the mysterious dark energy. 
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FIG. 5: The same as in Fig. |3] except for the Q m o — wo plane and the Q m o — w a plane. 
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